Continental Scale Hydrostratigraphy: Basin-Scale Testing of Alternative Data-Driven Approaches.

Integrated hydrological modeling is an effective method for understanding interactions between parts of the hydrologic cycle, quantifying water resources, and furthering knowledge of hydrologic processes. However, these models are dependent on robust and accurate datasets that physically represent spatial characteristics as model inputs. This study evaluates multiple data-driven approaches for estimating hydraulic conductivity and subsurface properties at the continental-scale, constructed from existing subsurface dataset components. Each subsurface configuration represents upper (unconfined) hydrogeology, lower (confined) hydrogeology, and the presence of a vertical flow barrier. Configurations are tested in two large-scale U.S. watersheds using an integrated model. Model results are compared to observed streamflow and steady state water table depth (WTD). We provide model results for a range of configurations and show that both WTD and surface water partitioning are important indicators of performance. We also show that geology data source, total subsurface depth, anisotropy, and inclusion of a vertical flow barrier are the most important considerations for subsurface configurations. While a range of configurations proved viable, we provide a recommended Selected National Configuration 1 km resolution subsurface dataset for use in distributed large-and continental-scale hydrologic modeling.

Continental Scale Hydrostratigraphy: Comparing Geologically Informed Data Products to Analytical Solutions.

This study synthesizes two different methods for estimating hydraulic conductivity (K) at large scales. We derive analytical approaches that estimate K and apply them to the contiguous United States. We then compare these analytical approaches to three-dimensional, national gridded K data products and three transmissivity (T) data products developed from publicly available sources. We evaluate these data products using multiple approaches: comparing their statistics qualitatively and quantitatively and with hydrologic model simulations. Some of these datasets were used as inputs for an integrated hydrologic model of the Upper Colorado River Basin and the comparison of the results with observations was used to further evaluate the K data products. Simulated average daily streamflow was compared to daily flow data from 10 USGS stream gages in the domain, and annually averaged simulated groundwater depths are compared to observations from nearly 2000 monitoring wells. We find streamflow predictions from analytically informed simulations to be similar in relative bias and Spearman's rho to the geologically informed simulations. R-squared values for groundwater depth predictions are close between the best performing analytically and geologically informed simulations at 0.68 and 0.70 respectively, with RMSE values under 10 m. We also show that the analytical approach derived by this study produces estimates of K that are similar in spatial distribution, standard deviation, mean value, and modeling performance to geologically-informed estimates. The results of this work are used to inform a follow-on study that tests additional data-driven approaches in multiple basins within the contiguous United States.

Water Table Depth Estimates over the Contiguous United States Using a Random Forest Model.

Water table depth (WTD) has a substantial impact on the connection between groundwater dynamics and land surface processes. Due to the scarcity of WTD observations, physically-based groundwater models are growing in their ability to map WTD at large scales; however, they are still challenged to represent simulated WTD compared to well observations. In this study, we develop a purely data-driven approach to estimating WTD at continental scale. We apply a random forest (RF) model to estimate WTD over most of the contiguous United States (CONUS) based on available WTD observations. The estimated WTD are in good agreement with well observations, with a Pearson correlation coefficient (r) of 0.96 (0.81 during testing), a Nash-Sutcliffe efficiency (NSE) of 0.93 (0.65 during testing), and a root mean square error (RMSE) of 6.87 m (15.31 m during testing). The location of each grid cell is rated as the most important feature in estimating WTD over most of the CONUS, which might be a surrogate for spatial information. In addition, the uncertainty of the RF model is quantified using quantile regression forests. High uncertainties are generally associated with locations having a shallow WTD. Our study demonstrates that the RF model can produce reasonable WTD estimates over most of the CONUS, providing an alternative to physics-based modeling for modeling large-scale freshwater resources. Since the CONUS covers many different hydrologic regimes, the RF model trained for the CONUS may be transferrable to other regions with a similar hydrologic regime and limited observations.

The evolution of dam induced river fragmentation in the United States.

It is established that dams decrease river connectivity; however, previous global scale studies of river fragmentation focused on a small subset of the largest dams. In the United States, mid-sized dams, which are too small for global databases, account for 96% of major anthropogenic structures and 48% of reservoir storage. We conduct a national evaluation of the evolution of anthropogenic river bifurcation over time that includes more than 50,000 nationally inventoried dams. Mid-sized dams account for 73% of anthropogenically created stream fragments nationally. They also contribute disproportionately to short fragments (less than 10 km), which is particularly troubling for aquatic habitats. Here we show that dam construction has essentially reversed natural fragmentation patterns in the United States. Prior to human development, smaller river fragments and less connected networks occurred in arid basins while today we show that humid basins are the most fragmented due to human structures.

ResOpsUS, a dataset of historical reservoir operations in the contiguous United States.

There are over 52,000 dams in the contiguous US ranging from 0.5 to 243 meters high that collectively hold 600,000 million cubic meters of water. These structures have dramatically affected the river dynamics of every major watershed in the country. While there are national datasets that document dam attributes, there is no national dataset of reservoir operations. Here we present a dataset of historical reservoir inflows, outflows and changes in storage for 679 major reservoirs across the US, called ResOpsUS. All of the data are provided at a daily temporal resolution. Temporal coverage varies by reservoir depending on construction date and digital data availability. Overall, the data spans from 1930 to 2020, although the best coverage is for the most recent years, particularly 1980 to 2020. The reservoirs included in our dataset cover more than half of the total storage of large reservoirs in the US (defined as reservoirs with storage greater 0.1 km3). We document the assembly process of this dataset as well as its contents. Historical operations are also compared to static reservoir attribute datasets for validation.

A hydrological simulation dataset of the Upper Colorado River Basin from 1983 to 2019.

This article presents a hydrological reconstruction of the Upper Colorado River Basin with an hourly temporal resolution, and 1-km spatial resolution from October 1982 to September 2019. The validated dataset includes a suite of hydrologic variables including streamflow, water table depth, snow water equivalent (SWE) and evapotranspiration (ET) simulated by an integrated hydrological model, ParFlow-CLM. The dataset was validated over the period with a combination of point observations and remotely sensed products. These datasets provide a long-term, natural-flow, simulation for one of the most over-allocated basins in the world.

Simulating Groundwater-Streamflow Connections in the Upper Colorado River Basin.

In mountain, snow driven catchments, snowmelt is supposed to be the primary contribution to river streamflows during spring. In these catchments the contribution of groundwater is not well documented because of the difficulty to monitor groundwater in such complex environment with deep aquifers. In this study we use an integrated hydrologic model to conduct numerical experiments that help quantify the effect of lateral groundwater flow on total annual and peak streamflow in predevelopment conditions. Our simulations focus on the Upper Colorado River Basin (UCRB; 2.8 × 105 km2 ) a well-documented mountain catchment for which both streamflow and water table measurements are available for several important sub-basins. For the simulated water year, our results suggest an increase in peak flow of up to 57% when lateral groundwater flow processes are included-an unexpected result for flood conditions generally assumed independent of groundwater. Additionally, inclusion of lateral groundwater flow moderately improved the model match to observations. The correlation coefficient for mean annual flows improved from 0.84 for the no lateral groundwater flow simulation to 0.98 for the lateral groundwater flow one. Spatially we see more pronounced differences between lateral and no lateral groundwater flow cases in areas of the domain with steeper topography. We also found distinct differences in the magnitude and spatial distribution of streamflow changes with and without lateral groundwater flow between Upper Colorado River Sub-basins. A sensitivity test that scaled hydraulic conductivity over two orders of magnitude was conducted for the lateral groundwater flow simulations. These results show that the impact of lateral groundwater flow is as large or larger than an order of magnitude change in hydraulic conductivity. While our results focus on the UCRB, we feel that these simulations have relevance to other headwaters systems worldwide.

Evapotranspiration depletes groundwater under warming over the contiguous United States.

A warmer climate increases evaporative demand. However, response to warming depends on water availability. Existing earth system models represent soil moisture but simplify groundwater connections, a primary control on soil moisture. Here we apply an integrated surface-groundwater hydrologic model to evaluate the sensitivity of shallow groundwater to warming across the majority of the US. We show that as warming shifts the balance between water supply and demand, shallow groundwater storage can buffer plant water stress; but only where shallow groundwater connections are present, and not indefinitely. As warming persists, storage can be depleted and connections lost. Similarly, in the arid western US warming does not result in significant groundwater changes because this area is already largely water limited. The direct response of shallow groundwater storage to warming demonstrates the strong and early effect that low to moderate warming may have on groundwater storage and evapotranspiration.

Simulating the sensitivity of evapotranspiration and streamflow to large-scale groundwater depletion.

Groundwater pumping has caused marked aquifer storage declines over the past century. In addition to threatening the viability of groundwater-dependent economic activities, storage losses reshape the hydrologic landscape, shifting groundwater surface water exchanges and surface water availability. A more comprehensive understanding of modern groundwater-depleted systems is needed as we strive for improved simulations and more efficient water resources management. Here, we begin to address this gap by evaluating the impact of 100 years of groundwater declines across the continental United States on simulated watershed behavior. Subsurface storage losses reverberate throughout hydrologic systems, decreasing streamflow and evapotranspiration. Evapotranspiration declines are focused in water-limited periods and shallow groundwater regions. Streamflow losses are widespread and intensify along drainage networks, often occurring far from the point of groundwater abstraction. Our integrated approach illustrates the sensitivity of land surface simulations to groundwater storage levels and a path toward evaluating these connections in large-scale models.

Connections between groundwater flow and transpiration partitioning.

Understanding freshwater fluxes at continental scales will help us better predict hydrologic response and manage our terrestrial water resources. The partitioning of evapotranspiration into bare soil evaporation and plant transpiration remains a key uncertainty in the terrestrial water balance. We used integrated hydrologic simulations that couple vegetation and land-energy processes with surface and subsurface hydrology to study transpiration partitioning at the continental scale. Both latent heat flux and partitioning are connected to water table depth, and including lateral groundwater flow in the model increases transpiration partitioning from 47 ± 13 to 62 ± 12%. This suggests that lateral groundwater flow, which is generally simplified or excluded in Earth system models, may provide a missing link for reconciling observations and global models of terrestrial water fluxes.